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Abstract 

The three photons from the dominant ortho-positronium decay and two vector 
mesons from the rj c exclusive decays are found to be in tripartite and high-dimensional 
entangled states, respectively. These two classes of entangled states possess the Hardy 
type nonlocality and allow a priori for quantum mechanics vs local realism test via 
Bell inequalities. The experimental realizations are shown to be feasible, and a con- 
crete scheme to fulfill the test in experiment via two-vector-meson entangled state is 
proposed. 
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Bell's theorem [T] states that no local hidden variable theory (LHVT) can reproduce all 
the possible results of quantum mechanics (QM). This can be testified via either inequalities, 
the Bell inequalities (Bis), satisfied by the expectation values of LHVT while violated by 
QM predictions, or logic contradictions remaining in LHVT vs QM. The original BI and 
its generalizations mainly concern about two-qubit systems and have been verified with 
various physical systems, see, e.g., [2j[3lll]. However, the analogous tests by virtue of high- 
dimensional and multipartite entangled states, especially in high energy experiment are very 
limited |3l [5] . 

In late 1980s, Greenberger, Home, and Zeilinger (GHZ) showed that certain three(or 
more)-partite entangled states may exhibit nonclassical effects more prominently than by 
bipartite state [6]. Soon afterwards, Hardy [7] invented a cute method to demonstrate the 
essential contradiction between QM and LHVT without employing the BI, in which the 
non-maximally entangled state works. A three-photon GHZ state was obtained in [8] by 
manipulating the parametric downconversion (PDC) photons and the experimental result 
deviates the local realistic expectation by more than eight standard deviations [9]. For the 
Hardy's type of violation, it was proved that the argument may take the form of Clauser and 
Home [10] (CH) inequality and then reduce the required detector efficiencies for loophole- 
free tests of Bell inequalities [TT1 [12] . Experimental measurement of the Hardy fraction was 
obtained in [13J. It is noteworthy that the existing entanglement tests are performed mainly 
in the regime of visible light. The non-maximally entangled kaon pairs were once suggested 
to realize the Hardy's type of criterion in experiment [HJ [15j [16] , however there has been no 
such experiment carried out yet to the best of our knowledge. 

In this Letter we propose to investigate the quantum non-local entanglement in two 
well-established high energy processes, in which the multipartite or/and high-dimensional 
entanglement is naturally realized, i.e., the three photons in the dominant ortho-positronium 
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decay and the two vector mesons in the r\ c exclusive decay. In [17], Acin, Latorre and 
Pascual already realized that the ortho-positronium decay can provide an entangled state 
of three space-separated photons. However, in this work we find that there is a Hardy type 
contradiction between QM and LHVT, and the CH type inequality would be violated by the 
quantum correlation of this state. 

In the middle of the last century, the polarization correlation of the two photons in para- 
positronium decay was used to determine its parity [TBI ED] . In quantum mechanics these 
two correlated photons are in the state of [21] 

l*> = ^{|*>l!/>-lv>l*>}, (i) 

where \x) and \y) stand for the photon polarization states in coordinate space. The ortho- 
positronium, the spin triplet state of electron and positron with J PC = 1 , decays over- 
whelmingly to three photons. In this case, the photon polarization information can be 
extracted from the decay amplitude, which is proportional to [22] 

V cc J2h( e 2 • e 3 - £ a • S 3 ) + <Ji(e 3 • 5 2 + e 2 • S 3 )} , (2) 

eye 

where e, denote the photon polarization vectors, ki represent their momenta, and c\ = fcj x e^. 
In the form of circular polarization, ef = t, L ± i5{ and the amplitude reads 

V<xJ2i4(e2 • 0+^(4-4)]- (3) 

eye 

Obviously, ef represent the left- and right-handed polarizations of photons, respectively. 
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Thus the three photon polarization state should be in the following form 



|37> 



1 



{\RRL) + \RLR) + \LRR) 



V6 



+ \LLR) + \LRL) + \RLL)} . 



(4) 



Here, L and R denote the left and right circularly polarized photons, respectively. 

As shown in |23j, under stochastic local operations and classical communication (SLOCC) 
pure three qubits entangled states have six inequivalent classes, and among them two are 
genuinely entangled, named as GHZ and W states. By virtue of the method proposed in 
|24j . one can easily confirm that the state (T4j) is SLOCC equivalent to GHZ state. On the 
other hand, from the point of entanglement measure, 3-tangle r defined in [25], the three- 
qubit entangled state (J4]) has r(|37)) = 1/3. Since r is none zero only in the case of GHZ 
or its equivalent states [23], we can also conclude that (jlj) is in SLOCC equivalence to the 
GHZ state. Therefore, we know that the three photons from the ortho-positronium decay 
constitute the simplest multipartite entangled system, the three-qubit state. 

To observe the non- locality of state (j3J), one needs also to measure the linear polarization 
of the state, e.g., the horizontal and vertical polarizations H and V. The transformation 
between circular and linear polarizations reads 



Under the linear polarization bases, the three-photon polarization wave function can be 
reexpressed as 





1 37) = -={3\HHH) + \HVV) + \VHV) + \VVH)} . 




(6) 
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To show the incompatibility of quantum mechanics with local realism(LR), we need an 
additional form of the state (jlj). Here, we express one photon in linear polarization, but left 
the other two in circular bases. In this case the wave function is 



|37> 



1 



{(\RR) + \LL) + 2\RL) + 2\LR))\H) 



Vl2 



i(\RR)-\LL))\V)}. 



(7) 



From Eqs.fll]), (jSJ) and (JZJ) one can easily get the following probabilities: 



^ 37 >(* = V,j = V) = -, P^iCj = C k \i = V) = l, 
P|3 7) (Ci = C k \ j = V) = 1 , P ]H (Ci = Cj = C k ) = . 
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Here, -P|3 7 }(« = V, j = V) denotes the probability of two of the photons being vertically 
polarized; P\ 3 ^(Cj = C k \i = V) denotes the conditional probability of photons j and k 
possessing the same circular polarization while the photon i being vertically polarized; and 
P\3y)(Ci = Cj = Ck) represents the probability of all three photons having the same circular 
polarization. 

The probabilities in (jSJ) can be embedded into CH type inequality [26], that is 



This inequality is a constraint imposed by the LHVT. While inputting the quantum me- 



P\^)(* = VJ = V)- P l3l) ^ = V, Cj + C k ) 



P\*r){Ci ^ C k ,j = V)- P^id = C 3 = C k ) < . 



(9) 
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chanics results on P|3 7 ), i.e., Eq. (jSJ), into ©, one readily finds that 

~-0-0-0<0, 
4 



MO; 



which is obviously a contradiction. Note that inequality (Q is valid only in the ideal case of 
perfect detection efficiencies, by which only a restricted class of Local Realistic Theories can 
be tested. 

To perform the experimental test on above inequality becomes realistic in recent years. 
Intense pencil beams of ortho-positronium (30-1.7 x 10 4 s _1 ) can be produced based on 
a low energy positron storage ring equipped with an electron cooling system [27J. This 
may accumulate about 10 12 three-photon events in ortho-positronium decay a year. Unlike 
in the case of Bell inequality where a group of measurements with relative angles in the 
measurement bases should be performed, here only two kinds of polarization measurements 
are necessary. However, we should point out that the measurement of photon polarization in 
positronium decay is a bit challenging since the photon energy lies in the regime of hard X-ray. 
The entangled photons from ortho-positronium decay are of hundreds of keV which is about 5 
orders greater than ordinary photons in visible region. In experiment, the X-ray polarimeter 
has made the measurements on linear and circular polarization components feasible in a 
wide range of energy. With ordinary inorganic crystal GaAs, people had already realized 
the measurement on polarizations of photon with energy span of 3-20 keV [28]. Recently, 
a novel kind of X-ray polarimeter based on hexagonal crystal, which can separate different 
polarization components simultaneously, is available [29J. Although its working energy is a 
bit lower than that of photon energy from positronium decay, further development on X-ray 
polarization analyzer is quite realistic in near future. 

As mentioned in the beginning, we find that a 3 x 3 entangled state in spin can be 
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produced in the process of r\ c exclusive decay into two light vector mesons, i.e., r\ c — > pp, (jxj), 
etc., as shown in Figure [TJ In the following part we briefly prove this, show its nature, and 
propose a scheme on how to carry out the test on LR by using of the two-vector-meson 
entangled state. Details will be presented elsewhere. 

According to the laws of angular momentum and parity conservation in strong interaction 
we know that those two vector mesons are in the state of total spin 1. From the Clebsch- 
Gordan decomposition, its zero component in z-direction is 

*|i,o> = ^(|l)|-l)-|-l)|l}),, (11) 

where the z-axis is chosen to be the moving direction of one vector meson in another meson's 
rest frame, as shown in Figure [2j Expressing the entangled state in bases perpendicular to 
z-axis, one has 

*H,o> = -^(|0> aj .|0> a -|0> a |0> a J. (12) 

Here, the subscript a indicates an arbitrary axis perpendicular to z, and a± means another 
base which is perpendicular to both z and a in terms of right-handed system. The plane 
fixed by axes a and a± is in fact the x-y plane, and a is the relative angle between x- and 
a-axes(or y- and a_i_-axes). Then one can immediately get four different probabilities within 
the x-y plane based on the QM, which are measurable quantities in experiment. They are 

P(J f3 = 0,J a = 0) = isin 2 (a-/3), (13) 

P(J^0,J 7 = 0) = Icos 2 (/3- 7 ), (14) 

P{ J x = 0, J a ^ 0) = icos 2 a, (15) 

P(J X = 0, J 7 = 0) = ^sin 2 7 , (16) 
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with 



P{Jp = 0,J a = 0) = | /3 (0U(0|^)| 1 , 0> | 2 , 
P(Jp ^ 0, J 7 = 0) = 

P(Jp = 1, J 7 = 0) + P(J^ = -1, J 7 = 0) , 
J Q = jjj. cos a + J y sin a . 

Here, P(r 1; r 2 ) means the probability under the conditions of 77 and r 2 , which are satisfied 
by the two partites of 1 and 2 of the entangled state, respectively. |0) Q stands for the 
eigenvector of angular momentum operator J a with eigenvalue J a = in the x-y plane. 
From the generalization of Hardy's argument to spin-s(s = systems [30J, we 

obtain the following LHVT constraint 

P(J x = 0,J 7 = 0) < P(J x = 0,J a ^0) 
+ P{ Jp ± 0, J 7 = 0) 

+ P( Jp = 0, J a = 0) . (17) 

Substituting the quantum mechanics results ffT3l) - f|T6|) into (fT7|) . one immediately obtains 
1 1 

-sin 2 7 < -[cos 2 a + cos 2 (/3 - 7) + sin 2 (a - ^)] . (18) 

This inequality is maximally violated by QM predictions while a = 3ir/8, f3 = it /A, 7 = 5ir/8, 
which gives a contradictive result of 2+ 8 v/2 < 6 ~g v/2 . 

In principle, in experiment one can determine the polarizations of ps or 0s through 
measuring the distributions of their two-body exclusive decay products, the its or Ks, and 
then find the difference between LR and QM predictions. Here, we find that there exists 
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another practical method to find out the difference in the experiment, that is to perform the 
measurement on the quantities in CH inequality event by event. Following we explain it in 
detail. 

Suppose the probabilities of a count being triggered by the decays of Vi, V 2 polarizing 
along Hi j2 being p(Hi, A) and q(H 2 , A) respectively, then we can readily get a CH inequality 

PHE2] 

P(H h n 2 ) - P(ni, n' 2 ) + P{ft 1 , H 2 ) + P{H[, H' 2 ) - P(ftj) - P(H 2 ) < . (19) 

Here, P(fii,n 2 ) = f p{n\, X)q(n 2 , A)p(A) dA, and hence P{n\) = f p(ni, A)p(A) dA with A 
being a set of hidden variables. 

In the quantum field theory, the differential decay width of r] c — > Vi(p,e*)V 2 (q,e'*) — > 
P(pi)P(p 2 )P(qi)P(q 2 ), as shown in Figure [21 takes the following form 



with 



oc k^k^i*)! 2 (20) 



KnxK^I*)! 2 = P(n 1; n 2 ), (21) 

where unit vectors Hi and H 2 are normalized projections of momenta pi and q\ respectively 
of the final pseudoscalar mesons in x-y plane, and 

l*) = ^(k,)K)-i%)K)) (22) 

with e X)V and e' x y being the transverse components of the polarization vectors e and e'. The 
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wave function of (1221) explicitly shows that the two vector mesons are entangled, the same 
as ( |T2i) . Note that our concerned process is similar to the B to two- vector-meson process, 
the B — > V1V2 [31], and hence the differential decay width (120]) can be analogously obtained. 
The details of how to derive (1201) and others will be presented elsewhere. 

Given that the four final pseudoscalars move with momenta p±, P2, qi, and the 
azimuthal angle tp between two decay planes of the entangled vector meson pair then equals 
to the angle between n\ and f?2, as shown in Figure EJ The magnitudes of P(ni,n2) in the 
CH inequality are therefore experimentally measurable, which is obviously the probability 
density, up to an overall normalization factor, from the definition of P(f?!,n 2 ). That is 

P^, n 2 ) = K ^ . (23) 

Here, k is a calculable normalization constant independent of specific theory for our aim, N 
is the total event number, and N(f) is the event number within azimuthal angle (p, the angle 
between n\ and In above expression, apart from the constant k the right-hand side is 
experimentally measurable, i.e., the differential decay width of r] c to four pseudoscalar mesons 
divided by its total width via intermediate vector meson 0s, and is hence also theoretically 
calculable in quantum theory. Inputting the theoretical and experimental results into (fT9l) . 
one may in principle find the incompatibility of quantum theory with LR. However, in 
practice, to perform the test of incompatibility the experiment efficiency should be taken 
into account. The general inequality efficiency and background levels was once discussed by 
Eberhard [TU [T2] , and for the wave function (T2"2"l) the violation of inequality ffTOl) yields the 
threshold efficiency r\ > 82.8% [33J. 

To carry out the test of Bell type inequality, the decay angles (ni, n2) should generally be 
chosen actively by experimenters, but this cannot be realized for mesons due to the passive 
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character of their decays |34j . Thus, here only a restricted class of LR can be tested, like 
in [35]. A genuine Bell test also requires the decay events of two vector mesons V\ and V2 
to be space-like separated. For the strongly decayed vector mesons (0, p, etc.), we cannot 
guarantee for each particular event of i] c — ► V V — ► {PP) {PP) that the decays of two vector 
mesons are space-like separated. The non-space-like decay events may induce the locality 
loophole in experiment. Fortunately the ratio of the space-like separated events can be 
determined by the magnitude of j3 [36]. And, for r\ c — > 00, it is easy to find that (3$ = 0.729 
[32] . which is greater than 0.59, the lower bound required to test the LR [35]. Note that 
the j3 has a negative effect on the threshold counter efficiency rj, and which then requires a 
higher detection efficiency in experiment. 

In summary, in this work we show that three photons in the decay of ortho-positronium, 
or other similar onium, form a tripartite entangled state, which is SLOCC equivalent to the 
GHZ state. Since in experiment the ortho-positronium decays to three-photon process can 
be well measured and huge data samples exist, new technology development on the X-ray 
polarimeter may enable a direct measurement on the hard photon polarization and hence 
find the difference between quantum theory and LR. Due to the detection efficiency and 
the unavoidable noise and losses, this cannot be considered as a loophole- free experiment. 
Nevertheless, the present method appears to be of broad interest because unlike the general 
production of GHZ state in optical cases the non-locality of the three photons in positronium 
decay arises in a dynamical process without postselection [37]. The investigation of three 
partite systems will show different features not only to that of bipartite qubit systems, but 
bipartite qutrit and in general bipartite qudit systems as well. 

We find that two vector mesons in r\ c exclusive decay automatically form a three-dimensional 
non-maximally entangled state. A concrete scheme to test the LR in experiment via the two- 
vector-meson entangled state is proposed. Since experimental measurements on exclusive two 
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body decay processes 



p -> mi, -> K + K-(K° L K° S ) (24) 

are well-established, and give large branching fractions of p — > 7r7r ~ 100% and — > 
A+A~(A£Ag) ~ 49.2 ± 0.6% (34.0 ± 0.5%) [38J, we believe that in current running ex- 
periments the CH type inequality can be readily measured. It should be mentioned that the 
non-space-like decay events may induce the so-called locality loophole in experiment, which 
hinders the proposed tests to refute the LR definitely. Nevertheless, the experimental real- 
ization of the proposals in this work may extend the test on local realism into high energy 
regime with multipartite and high dimension, which will give us a more explicit conclusion 
in comparison with that from the bipartite qubit results. 
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Figure Captions 



Figure 1: The schematic Feynman diagram of processes r] c — > VV. Here V stands for vector 
meson <p or p. 

Figure 2: r\ c — > Vi(p, e*)V 2 (g, e'*) — > [P(pi)P(p 2 )][-P(?i)-P(?2)] decay kinematics in the rest 
frame of Vi- Here, Vi, ^2 stand for vector meson pair of 0s or ps, -P(pj) and P(qi) for the 
pseudoscalar mesons in V\ and V 2 decays, the Ks or 7rs . 
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Figure 1: 
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